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Trade  name  is  used  in  this  publication  solely  for 
the  purpose  of  providing  specific  information. 
Mention  of  a  trade  name  does  not  constitute  a 
guarantee  or  warranty  of  the  product  by  the 
Federal  Government  or  an  endorsement  by  the 
Government  over  other  products  not  mentioned. 


COMPUTER  REDUCTION  OF  PRECIPITATION  DATA 

George  H.    Comer  and  W.    Russell  Hamon 


Abstract 

A  program  has  been  developed  for  the  IBM   1620  computer  for  reduction  of 
digitized  precipitation  data  into  useful  quantities  for  analysis  and  reporting. 
The  input  format  used  with  the  program  is  designed  to  handle  manually  tran- 
scribed data  or  data  punched  by  an  automatic  chart  reader  from  Weather 
Bureau-type  recording  rain  gage  charts.    For  maximum  flexibility,    specific 
sections  of  the  program  may  be  bypassed.    Each  type  of  output  has  a  code 
unique  to  that  group. 

Flow  diagrams,    program  listings,    and  complete  instructions  for  using  the 
program  are  included. 


Introduction 


The  reduction  of  precipitation  data  from  the  analog  form,    as  obtained 
from  the  Weather  Bureau-type  recording  rain  gage,    to  quantities  needed  for 
analysis  and  reporting  is  a  tedious  and  time-consuming  task.    Such  reduced 
data  are  essential  in  the  study  of  runoff  rates,    water  yields,    sedimentation, 
and  other  hydrologic  research.    In  a  well-organized  and  comprehensive  pro- 
gram of  such  hydrologic  investigations,    large  volumes  of  precipitation  data 
must  be  processed;  and  the  electronic  computer  has  great  utility  in  such  a 
task. 

This  publication  describes  a  computer  program  developed  to  enable  the 
IBM     1620  computer  to  perform  all  the  calculations  required  in  reducing  the 
precipitation  data.    The  program  is  basically  the  same  as  that  developed  by 
Comer.,2    with  some  additions  and  refinements  incorporated  in  the  new  pro- 
gram.   Although  the  program  described  herein  was  written  for  the  IBM   1620 
computer,    the  methods  and  principles  used  could  be  employed  to  develop 
similar  programs  for  other  computers. 

Agricultural  engineer  Northeast  Branch,  Danville,  Vt. ;  and  research  agricultural  engineer,  Southern  Branch, 
Oxford,   Miss. ,   Soil  and  Water  Conservation  Research  Division,  Agricultural  Research  Service,  USDA,  respectively. 

Comer,    George  H.    Computer  analysis  of  rainfall.  Unpublished  M.  S.  thesis,  Newman  Library,  Virginia 
Polytechnic  Institute,  Blacksburg,  Va. ,  1961. 
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Data  Processing  Requirements 

The  source  of  the  raw  data  and  the  input  and  output  formats  were  basic 
considerations  in  writing  the  computer  program.    The  input  format  was  de- 
signed to  provide  an  easy  means  of  transcribing  the  analog  or  chart  data  by 
both  manual  and  automatic  means  into  a  form  usable  by  the  computer.    The 
output  format  was  arranged  to  contain  all  the  desired  reduced  precipitation 
data  in  a  form  easily  understood  for  direct  use  and  appropriately  identified 
for  use  as  input  data  to  computer  analysis  programs. 

Another  important  consideration  was  program  flexibility.    The  largest 
single  expense  in  reducing  the  precipitation  data  is  the  cost  of  computer  time. 
Any  time  used  in  computing  data  not  required  is  an  unnecessary  expense.    To 
eliminate  such  expenses  and  the  need  for  several  programs,    the  program  was 
written  with  provision  to  bypass  specific  sections. 

The  following  were  used  as  the  basis  for  the  program: 

1.  Compute  the  intensity  and  amount  of  precipitation  for  each  time 
interval  during  a  storm  when  the  intensity  may  be  assumed  as  con- 
stant. 

2.  Compute  running  totals  of  precipitation  for  each  storm,    day,    and 
month. 

3.  Compute  the  time  between  storms,   with  storms  determined  accord- 
ing to  item  7  below. 

4.  Compute  for  each  storm  of  0.25  inch  or  more  of  precipitation  the 
maximum  intensities  of  precipitation  during  the  storm  for  periods 
of  2,    5,    10,    15,    20,    and  30  minutes  and  1,    2,    4,    6,    12,    and  24 
hours. 

5.  Compute  the  amounts  of  precipitation  for  any  three  of  the  regular 
time  periods  of  1,    2,    3,    5,    6,    10,    15,    30,    or  60  minutes. 

6.  Compute  the  total  time  of  actual  precipitation  for  each  month  and 
for  each  year. 

7.  Separate  storms  by  periods  of  1  hour  or  more  during  which  no  pre- 
cipitation occurred,  with  provision  to  change  this  requirement  from 
storm  to  storm  by  input  coding. 

8.  Eliminate  the  necessity  to  start  each  period  of  precipitation  at  zero 
chart  depth  reading  by  providing  for  direct  reading  of  rain  gage  re- 
corder charts. 

9.  Provide  for  transcribing  the  input  data  by  manual  means  or  by  au- 
tomatic chart  reader. 

10.     Code  and  number  the  output  data  to  facilitate  automatic  handling  of 
the  data. 


The  Computer  Program 


The  Symbolic  Programing  System  (SPS)  was  chosen  as  the  language  to  be 
used  in  writing  the  precipitation  program  because  it  gives  a  flexible  program 
that  is  easily  altered.    It  is  also  the  easiest  method  of  writing  in  fixed  point 
arithmetic. 

A  complete  listing  of  the  program  source  deck  is  included  as  appendix  1; 
the  compressed  precipitation  program  deck  is  listed  as  appendix  2;  and  the 
flow  diagram  is  shown  as  appendix  3. 

A  standard  IBM   1620  computer  with  20,  000-digit  storage  capacity  and  au- 
tomatic division  feature  can  process  a  storm  with  a  maximum  of  1,  260  lines 
of  input  data  covering  a  maximum  duration  of  166  hours  and  38  minutes. 

The  input  format  used  with  this  program  is  designed  to  handle  manually 
transcribed  data  or  data  jpunched  by  an  automatic  chart  reader  with  a  freely 
variable  X  origin  offset  and  scale  control  ranging  from   -999  to  +999.    Thirty 
columns  (thirty-one  with  12 -hour  charts  as  punched  with  the  automatic  chart 
reader)  of  input  data  are  used.    Figure   1  shows  a  sample  form  for  manual 
listing  of  the  input  data.    The  keypunch  operator  punches  these  data  in  input 
cards. 
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FIGURE  1.  --Sample  form  for  manual  listing  of  input  data;  K  column  designates  separation  of  storms; 
M  column  designates  method  used  to  transcribe  input  data. 


The  input  card  formats  with  sample  data  punched  are  shown  in  figure  2. 
Columns  1  to  5  contain  the  watershed  number.    The  first  two  columns  represent 
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<FIGURE  2.  --Input  data  cards  with  same  data:  A,  Punched  from  manually  transcribed  data;  B,  punched 
by  a  cliart  reader  from  a  24-hour  chart;  and  C,  punched  by  a  chart  reader  from  a  12-hour  chart. 


"location"  identifiers,    corresponding  to  the  page  Nos.  of  the  USDA-ARS  Hy- 
drologic  Data  Releases,    with  columns  3  to  5  denoting  the  watershed  in  which 
the  rain  gage  is  located  or  situated  nearest.    Columns  6  to  8  contain  the  rain 
gage  No.  ,   and  columns  9  to   12  contain  the  card  number  which  may  be  ignored 
but,    when  punched,    affords  an  easy  means  of  restoring  the  input  data  to  its 
proper  order  when  disturbed. 

The  date  (month,    day,    and  year)  of  the  reading  is  entered  in  columns   13  to 
18  and  the  time  of  day  in  columns   19  to  ZZ.    If  the  data  are  manually  tran- 
scribed and  punched,    the  time  columns  contain  the  time  of  day  in  military 
time.    In  case  the  automatic  chart  reader  is  employed  with  the  Z4-hour  re- 
cording rain  gage  chart  as  shown  in  figure  3,    the  setup  below  is  used. 

Z4-hour  Rain  gage  Charts  (fig.    Z,    B) 

Time  Columns 


1 19  20  21  22 

8:00  a.  m -  9  6  0 

1Z:00  noon 7  Z  0 

l:00p.m -  6  6  0 

1Z:00  midnight 0  0  0 

12:01  a.m b  0  0  1 

7:59a. m b  4  7  9 


The  symbol  "b"  denotes  blank  column. 

For  the   lZ-hour  recording  rain  gage  chart  where  the  time  of  chart  begin- 
ning is  noon  and  midnight,    the  setup  below  is  used. 

lZ-hour  Rain  gage  Charts  (fig.    Z,    C) 
Time  Columns 
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The  chart  line  depth  in  inches  is  contained  in  columns  23  to  26.    This  depth 
is  taken  at  the  time  recorded  in  columns    19  to  22.    The  decimal  shown  in  fig- 
ure  1,    under  chart  line  depth,    is  understood  and  not  punched. 

The  K  columns  (columns  27  and  28)  of  the  card  are  used  to  designate  the 
separation  of  storms  according  to  the  following  coding  system: 

07--The  computer  will  consider  this  reading  as  the  first  reading  of 

another  storm  regardless  of  the  time  interval  between  reading  or 
the  amount  of  precipitation  occurring  since  the  previous  reading. 
This  code  causes  the  computer  to  consider  the  previous  reading  as 
the  end  of  a  storm. 

08--The  computer  separates  the  events  as  programmed.    (Any  two  periods 
of  precipitation  separated  by  1  hour  or  more  with  no  additional 
precipitation  between  are  considered  as  separate  storms.  ) 

09--The  computer  will  consider  this  reading  as  part  of  the  same  storm 
regardless  of  the  time  interval  or  amount  of  precipitation  occur- 
ring since  the  previous  reading. 

The  M  columns  (columns  29  and  30)  are  used  to  designate  the  method  used 
to  transcribe  the  input  data.  An  08  code  indicates  input  from  a  chart  reader;  a 
09  code  indicates  the  data  were  manually  transcribed  and  punched. 

The  output  data  consist  of  five  types:  (1)  Output  for  each  input  reading, 
(2)  extra  output  giving  amounts  of  precipitation  for  regular  time  periods,    (3) 
maximum  storm  intensities,    (4)  monthly  summaries,    and  (5)  yearly  summa- 
ries.   Each  type  of  output  has  a  code  unique  to  that  group. 

The  output  data  are  punched  in  logical  order  and  each  op.tput  card  is  num- 
bered to  facilitate  automatic  sorting  and  other  automatic  data  handling  pro- 
cedures. 

The  program  may  be  considered  as  a  combination  of  three  separate  parts. 
Provision  is  made  for  using  the  following  combinations  of  the  three  parts: 

1.  Part  1  only.  3.    Parts   1  and  3. 

2.  Parts   1  and  2.  4.    Parts   1,    2,    and  3. 

Part  1  of  Program 

Part   1  of  the  program  computes  the  intensity  and  amount  of  precipitation 
for  each  time  interval  and  accumulates  the  total  precipitation  for  each  storm, 
day,    and  month.    It  also  computes  monthly  and  yearly  summaries  and  the  num- 
ber of  hours  of  no  rain  between  storms. 

The  output  of  part  1  contains  the  following  data: 

1.  Watershed  No.    and  rain  gage  No.  ,    columns   1  to  8. 

2.  Date  (month,    day,    and  year),    columns  9  to   14. 
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3.  Time  of  day  (military  time),    columns   15  to  20. 

4.  Time  interval  (hours  and  minutes),    columns  21  to  28. 

5.  Intensity  for  the  time  interval  in  inches  per  hour  accurate  to  four  deci- 
mal places,    columns  29  to  34. 

6.  Amount  of  precipitation  for  the  time  interval  in  inches,    accurate  to  two 
decimal  places,    columns  35  to  40. 

7.  Accumulated  amounts  of  precipitation  for  storms,    days,    and  months  in 
inches  accurate  to  two  decimal  places,    columns  41  to  58. 

8.  A  code  (0--numeric  zero)  as   "blank"  from  IBM  407  accounting  machine, 
column  59. 

9.  The  sequential  card  number,    columns  75  to  78. 

Part   1  of  the  program  leaves   17  columns  of  the  standard  80-column  card 
unused.    Fifteen  of  these  columns  (60  to  74)  are  used  for  the  output  from  part 
3  when  it  is  used.    The  two  remaining  columns  (79  and  80)  are  not  used. 

At  the  end  of  each  month  a  summary  card  is  punched.    It  contains  the  data 
indicated  below: 

1.  Watershed  No.    and  rain  gage  No.  ,    columns   1  to  8. 

2.  Month  and  year,    columns  9,    10,    13,    and  14. 

3.  Time  of  actual  precipitation  during  the  month  (hours  and  fractions  of 
hours),    columns  21  to  28. 

4.  The  total  amount  of  precipitation  for  the  month  (inches),    columns   53 
to  58. 

5.  A  code  (6)  column  59. 

6.  The  sequential  card  number,    columns  75  to  78. 

A  similar  summary  card  is  punched  at  the  end  of  each  group  of  data  proc- 
essed (yearly  summary  card).    This  card  contains  the  same  type  of  informa- 
tion.   The  time  duration  of  precipitation  and  the  amount  of  precipitation  are 
for  the  period  covered  by  the  data  being  processed.    The   "year"  only  is  listed 
in  the  date  columns. 

The  first  output  card  punched  at  the  beginning  of  each  storm  contains  the 
number  of  hours  since  the  end  of  the  previous  storm  or  the  hour  of  no  rain 
between  storms. 


Part  2  of  Program 

Part  2  of  the  program  determines  maximum  precipitation  intensities  for 
each  storm  of  0.  25  inch  or  more  of  precipitation.    These  storms  are  called 
storm  events.    The  following  information  is  punched  in  these  output  cards: 

1.  Watershed  No.    and  rain  gage  No.  ,    columns   1  to  8. 

2.  Date  of  storm  ending  (month,    day,    and  year),    columns  9  to   14. 

3.  The  maximum  intensities  for  2,    5,    10,    15,    20,    and  30  minutes,    all  in 
inches  per  hour  accurate  to  four  decimal  places,    columns  2  1  to  58. 

4.  A  code  (9),    column  59. 

5.  The  sequential  card  number,    columns  75  to  78. 

6.  The  storm  event  number,    columns  79  and  80. 

A  second  card  is  punched  containing  the  same  data  except  the  maximum 
intensities  are  for   1,    2,    4,    6,    12,    and  24  hours,    columns  2 1  to  58. 

Part  3  of  Program 

Part  3  of  the  program  computes  precipitation  amounts  for  three  regular 
time  periods.  The  program  is  written  so  that  each  time  period  will  have  an 
ending  time  every  hour,    clock  time.    For  example: 

Consider  a  rainfall  that  begins  at   1:38  a.m.    and  the  three  regular  time 
intervals  being  used  are   15,    30,    and  60  minutes.    The  time  intervals  are  then 
as  follows: 

15-minute  intervals -- 1:  30- 1:45  30-minute  intervals -- 1:30 -2:00 

1:45-2:00  2:00-2:30 

2:00-2:15  2:30-3:00 

60-minute  intervals --1:00-2:00 

2:00-3:00 
3:00-4:00 

The  output  for  this  part  of  the  program  occupies  15  of  the  card  columns 
not  used  by  part  1  (columns  60  to  74).  The  data  contained  in  these  columns 
are  as  follows: 

1.  The  amount  for  the  first  regular  time  period  in  inches  accurate  to 
three  decimal  places,    columns  60  to  63. 

2.  The  amount  for  the  second  regular  time  period  in  inches  accurate  to 
three  decimal  places,    columns  64  to  67. 

3.  The  amount  for  the  third  regular  time  period  in  inches  accurate  to  two 
decimal  places,    columns  68  to  70. 


4.    The  time  of  ending  of  period  (military  time),    columns  7  1  to  74. 

If  no  regular  time  period  is  completed  during  the  time  interval  between 
readings,    these  columns  are  not  used.    When  more  than  one  of  any  regular 
time  period  can  be  accumulated  between  chart  readings,    extra  cards  are 
punched  containing: 

1.  Watershed  No.   and  rain  gage  No.  ,  columns   1  to  8. 

2.  Date  (month,    day,    and  year),    columns  9  to   14. 

3.  Time  of  day  of  period  ending  (military  time),    columns   15  to  20. 

4.  A  code  (5),    column  59. 

5.  The  precipitation  amounts  for  the  three  regular  time  periods,    columns 
60  to  70. 

6.  The  time  of  ending  of  period  (military  time),    columns  71  to  74. 

7.  The  sequential  card  number,    columns  75  to  78. 

At  the  end  of  each  storm  an  extra  card  of  the   same  form  is  punched  giving 
the  accumulation  in  each  regular  time  period  from  the  previous  period  until 
the  storm  ending. 

The  program  as  listed  in  appendix  2  will  compute  the  precipitation  amounts 
for  the  regular  time  periods  of  15,    30,    and  60  minutes.    To  compute  the 
amounts  for  the  regular  time  periods  of  5,    10,    and   15  minutes,    the  card  num- 
ber  "10"  of  the  compressed  program  deck  listing  (appendix  2)  must  be  re- 
placed by  a  card  of  the  form  shown  in  figure  4,    A_.    Similarly,    computations 
for  time  periods  of  2,    5,    10  minutes  may  be  accomplished  by  substituting  a 
card  of  the  form  shown  in  figure  4,    B. 

A  sample  of  the  output  listed  by  an  IBM  407  accounting  machine  is  shown 
in  figure   5.    To  list  this  data  the  panel  board  for  the  407  must  be  wired  for  the 
data.    The  wiring  diagram  for  the  panel  board  is  included  as  appendix  4. 

Header  Card 

The  header  card  (fig.    6)  is  used  to  set  up  the  output  area. 

When  part  of  a  year  of  data  is  processed,    a  different  header  card  should 
be  used  in  computing  the  remaining  data  for  the  year.    The  header  card  for 
this  case  would  be  punched  as  follows: 

1.  Duplicate  the  output  card  corresponding  to  the  last  input  card  except 
the  card  No. 

2.  Make  the  card  No.   the  same  as  the  card  No.    of  the  last  output  card. 

3.  Use  this  card  as  the  header  to  process  the  remaining  data  in  that  year. 
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Figure  4.  --Form  of  program  card  required  to  compute  precipitation  amounts  for  regular  time  periods 

of:  A,  5,   10,  and  15  minutes;  B,  2,  5,  and  i0  minutes. 


Trailer  Card 

The  trailer  card  will  be  punched  in  either  of  two  ways,    depending  on  the 
last  input  card. 

If  the  last  input  card  is  the  last  reading  of  a  storm,    the  trailer  card  will 
be  blank  with  the  exception  of  "9"  punches  in  columns  79  and  80  (type   1  trailer 
card).   When  the  last  input  card  is  not  the  end  of  a  storm,    the  trailer  card  is 
blank  except  column  80  contains  a  "9"  punch  (type  2  trailer  card). 

First  Reading  in  the  Year 

The  first  reading  in  each  group  of  data  being  processed  must  have  a  chart 
line  depth  reading  of  zero. 
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Figure  5.  --Sample  of  output  listed  by  IBM  407  accounting  machine. 


Last  Reading  in  the  Year 

When  there  is  no  storm  in  progress  at  the  end  of  a  year  and  the  number  of 
hours  since  the  last  precipitation  is  desired,  the  last  card  should  be  for  mid- 
night of  December  31  with  a  chart  line  depth  reading  of  zero.  In  this  case  use 
type   1  trailer  card. 

In  case  a  storm  crosses  a  year  the  procedure  below  may  be  used. 

1.    Take  the  last  reading  in  the  year  at  midnight  December  31.    Follow  this 
by  a  type  2  trailer  card  and  two  blank  cards. 


2.    Process  this  year  of  data  according  to  "Instructions  for  Operation"   1  to  9, 
page   14. 
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Figure  6.  — Sample  header  card  used  to  set  up  output  area. 

3.  Immediately  after  processing  the  above  data,    process  the  data  in  the 
next  year  (at  least  until  the  end  of  the  storm  crossing  into  the  following 
year). 

4.  The  header  card  for  the   "new  year"  would  be  punched  as  shown  in  fig- 
ure 7,    except  the  SAM  columns  that  would  contain  the  storm  accumu- 
lation as  of  midnight  December  31  of  the   "old  year.  "  The  first  reading 
in  the   "new  year"  must  be  taken  as  if  the  chart  line  depth  were  zero  at 
midnight  December  31  of  the  "old  year.  " 

An  alternate  procedure  would  treat  the  storm  crossing  the  year  as  two 
storms,    i.e.,    one  storm  ending  at  midnight,    December  31,    and  the  other  be- 
ginning at  this  time.    To  process  the  data  using  this  procedure,    use  the  steps 
presented  below. 

1.  Take  the  last  reading  in  the  year  at  midnight  December  31.    Follow  this 
by  a  type   1  trailer  card  and  two  blank  cards. 

2.  Process  this  year  of  data  according  to  "Instructions  for  Operation"   1 
to  9. 

3.  The  last  part  of  the  storm  would  be  processed  with  the  next  year  of 
data. 

4.  The  first  reading  in  the  "new  year"  should  be  taken  at  1  minute  after 
midnight  with  the  chart  line  depth  equal  to  zero  and  other  readings  in 
the  storm  similarly  adjusted.    The  remaining  data  in  the  year  would  be 
handled  in  the  normal  manner. 

Speed  of  Operation 

The  speed  of  operation  will  depend  on  the  parts  used.    Part  1  alone  oper- 
ates at  punch  speed,    or   125  readings  per  minute.    Parts   1  and  3,    when  used 
together,    operate  at  punch  speed.    However,    the  extra  output  cards  for  part  3 
reduce  the  number  of  readings  processed  per  minute. 
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In  processing  36  station  years  of  data  from  Danville,    Vt.  ,    when  parts    1 
and  2  were  used,    an  average  of  25  minutes  was  required  to  process  a  year  of 
data.    The  input  decks  averaged  800  cards  and  covered  an  average  of  42  storm 
events  per  year.    For  all  three  parts  used  at  Oxford,    Miss.  ,    an  input  deck  of 
850  cards  required  40  minutes.    The  output  contained  2,  000  cards  covering  56 
storm  events. 


Instructions  for  Operation 

1.  Clear  memory. 

2.  Place  the  compressed  program  deck  followed  by  two  blank  cards  in  the 
reader  hopper  and  press   "LOAD.  " 

3.  After  the  last  card  is  read  in,    depress   "NON  PROCESS  RUNOUT"  until 
the  two  blank  cards  come  out. 

4.  Set  the  console  switches  as  follows: 

Parity  -   "STOP" 

1/0  -  "STOP" 

O-Flow   -   "STOP" 

Program  switch  1   -  "OFF.  "  To  omit  the  maximum  intensity  deter- 
mination (part  2),    place  switch  in  "ON"  position. 

Program  switch  2  -  "OFF.  "  To  omit  the  determination  of  precipita- 
tion amounts  for  regular  time  periods  (part  3),  place  switch  in  "ON" 
position. 

Program  switch  3--Not  interrogated. 

Program  switch  4--Not  interrogated. 

5.  Place  blank  cards  in  the  punch  hopper. 

6.  Place  the  input  data  in  the  reader  hopper  in  the  following  order: 

a.  Header  card 

b.  Data  cards 

c.  Trailer  card 

d.  Two  blank  cards 

7.  Press   "READER  START"  and  "PUNCH  START.  " 

8.  Press   "START.  " 

9.  After  the  last  card  has  been  processed  the  computer  will  stop  with  the 
"READER  NO  FEED"  light  on.  Press  both  "NON  PROCESS  RUNOUT" 
switches  and  remove  all  cards. 
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10.    If  you  have  another  year  of  data  to  run,    repeat  steps  4  to  10  and  omit 
step  8. 


Summary 

The  tedious  and  repetitious  calculations  required  in  reducing  raw  precipi- 
tation data  into  quantities  needed  for  analysis  and  reporting  can  be  performed 
by  an  IBM   1620  computer  using  the  program  as  presented.    Input  data  for  the 
program  may  be  punched  either  manually  or  by  an  automatic  chart  reader. 

Precise  calculations  are  performed  by  the  computer  in  determining  values 
for  three  separate  output  sections:  namely,    (1)  precipitation  intensities  and 
amounts  of  precipitation  by  intervals,    storm,    day,    month,    and  year;  (2)  max- 
imum intensities  for  various  time  periods;  and  (3)  precipitation  amounts  for 
regular  time  periods.    The  last  two  sections  may  be  bypassed  as  desired. 

The  output  data  are  punched  in  cards  of  a  form  designed  for  re -use  as 
"input"  data  to  analysis  and  data  handling  programs.    Each  card  is  numbered 
and  each  type  of  output  has  a  code  unique  to  that  group,    so  that  all  sorting, 
separation,    and  handling  of  the  data  may  be  done  by  machine. 
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APPENDIX  1 


SPS  SOURCE  DECK  LISTING  FOR  IBM  1620  COMPUTER 


D0RG402 
TTFC  DC   2,0 


PMTF   B    AROUND-12,, 
ORG    BNC1RETURN,,0 


0 


N1 
N2 


DS 
DS 


IN3 
IN4 

DS 
DS 

4 
2 

IN5 
IN6 

DS 
DS 

2 
2 

IN7 

IN8 

DS 
DS 

4 
4 

IN9 
INI  1 

DS 
DS 

2 
2 

IN12 

DS 
DS 

4 
18 

0UT1 

DS 
DS 

28 
5 

0UT2 
0UT3 

DS 
DS 

3 

2 

0UT4 
0UT5 

DS 
DS 

2 
2 

0UT6 

0UT7 

DS 
DS 
DS 
DS 

6 
8 

OUTB 
0UT9 

6 

0UT10 
0UT11 

D$ 
DS 

6 
6 

0UTT2 
0UT13 

DS 
DS 

7 
4 

OUT  14 

0UT15 

DS 
DS 

4 
3 

0UT16 
0UT17 

DS 
DS 

4 
4 

0UT18 
EOUTl 

DS 
DS 

2 

5 

E0UT2 
E0UT3 

DS 
DS 

3 
2 

E0UT4  DS 
E0UT5  DS 


0UT19 
0UT20 


0UT21 
_gUT_22 
OUT23 
0UT24 


DS 

DS 

DS 

DNB 

DS 

DNB 


14 

6 

8" 

7 
15 


0UT25  DS  W 

0UT26  DS  2 

CODES  DC  ZT76789~ 

N01    DC  4A0 

NO 2   DC  4,0 

NOMO   DC  2,0 


E0UT6  DS 
E0UT7  DNB 

4 
40 

E0UT8  DC 
E0UT9  DS 

5,50000 
4 

E0UT10DS 
E0UT11DS 

3 

4 

E0UT12DS 
E0UT13DS 

4 
2 

JAN 

fEB_ 

MAR 

APR 

MAY 

JUNE 


ZERO   DC   10,0 
-X5- fl£__A,H.._ 

X6    DC   4,0 
Ml D.C_._a,Q._ 

DE3   DC   9,0 
-K2 DC   9,0 


TIME   DC 

YIIMELDJI. 

YAMT  DC 

JTFJJ D£_ 


TIFH 
XE 


DC 


8,0 
6,0 
11,0 


A 

B 

AMT1 

AMT2 

AMT3 

TCK1 


DC 
DC 
DC 
DC 
DC 
_D_C_ 


4,0 

6,0 

6,0 
2.0 


TCK2   DC 

...ICK3— -D.C. 

T.I  1    DC 

T J  2 DC 

TI3   DC 
DIGIT  DC 


2,0 
JUQ  . 

2,15 
-2,30 


2,60 
2JL 
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APPFNiniX  1  mNTINUED 


REK1 


IPM   DC 

_CJ4Ifi_D£- 
DC 

.DC 
DC 
DC 
DC 
DC 


6,0 


4,0 

-4.0- 

3,0 

4f@ 


B3 

TIA   DC 

.AM¥__ Q£ 

ROAH   DC 
TIAH DX_ 


7,0 

4,0 

7,0 
-AJL 


LD   92,0UT6-2 
DM  91.60,10 


START 


TF   TFH,97 

A   TFH-7.0UT6-4 


RNCDlTTT^f 
SF   INI -4 


SF 
SF 


N2-2 
N4-1 


SF 
SF 


HOLD 

CNT_ 

0UT30 

-01IT3.1- 
OUT32 


0UT34 

OUT  3  5 

.QU.T36 

OUT37 

.OUT38, 


DC   6,0 
J)SA_SIQRE=24 
DC   5,0 

DC.  _5^D 

DC   2,0 

JD_C Z^Q 


DC   2,0 
DNB_a._._ 


DC 
HC 
DC 


OUT39 
QUT40 

DUT41 

0UT42 

-1NT1- 


DC 
DC 
DNB 
DC 
DC 
_QC_ 


6,0 
6,0 
6,0 


SF 
SF 


N5-1 
N6-1 


SF   IN11-1 

BD   *+24,0UT1-5 


B 
BD 


*+72 
*+24,0UT1-6 


B   MOUT 
BNC2F0UT 


BNC1SEARCH 
B    MOUT 


C 
BNN 


l"NH7'X5 
NEXT 


TF 
NEXT   C 


DC 
DC 
DC 
DC. 
DC 
_D_C_ 


6,0 
7,0 
15 

2,0 


DC 
_D_C_ 


<*,5 
4,10 
4,15 
A..20. 
4,30 

AM. 


X5, IN8 
IN4,0UT3 


BE   NEXT1 

C OUT3.ZERO-8 


BE   DELTA 
BNC1DELTA 
BNC2DELTA 
MOUT   TF   0UT19,0UT5 


DC 

^DC 

DC 

DC 


4,120 
k+2kQ 


4,360 
4,720 
4,1440 
4,9999 


TF      0UT19-2,ZER0-8 

IF  „0UT19-4,0UT3 
TF  6UT19-6,bUT2 
TF      0UT19-9.0UT1 


TF   0UT21,TIME 

T  F      0UT23-1.0UT12 - 1 


TD 
CF 


RECORDDC      6,0 
DC      6.^0 


DC 

DC 
DC 


~67o 

6.0 


b/0 

1&. 


AM 
TF 


OUT23,CODES-3 
0UT23  _ 

0UT17, 1 ,8 
OUT25.0UT17 


WNCD0UT19-13 
A        YTIME.TIME 


A 
_TF 

TF 
TF 


YAMT,0UT12-1 
T IME.ZERO-2 


OUT12,ZERO-3 
0UT11.ZER0-4 


FIRST  RNCD0UT1-4 
CF   E0UT8-4 


BD   END,0UT1-5 

.  ,   BNC1RETURN 

BNC2RETURN-24 
DELTA  C   OUT6.ZERO-6 
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APPENDIX-  1  CONTINUED 


BNE 
CM 

*+48 
0UT6-4,0,10 

BNE 
B 

*+24 
NEXT2 

TFM 
S 

DE3-7,24,10 
DE3.TFH 

NEXT2 

TF 
LD 

TFH,ZER0-1 
98, IN6 

DM 
C 

98,40,10 
99.ZER0-8 

BE 

TFM 

NEXT3 
FEB, 28, 10 

NEXT3 

B 
TFM 

NEXT4 
FEB, 29, 10 

NEXT4 

MM 
A 

0UT3,2,9 
NEXT5+11.99 

MM 
SM 

IN4,2,9 
99f2 

NEXT5 

AM 
A 

99,N0M0 
M1.N0M0.7 

C 
BNE 

NEXT5+11,99 
*+36 

TFM 
B 

NEXT5+11,N0M0 
NEXT6 

AM 
B 

NEXT5+1 1 , 2 

NEXT5 

NEXT6 
NEXTl 

SM 
CM 

Ml, 1,9 
IN11,9,10 

BNE 
SF 

READER 
IN7-3 

SF 
LD 

IN7-1 
92. IN7 

DM 
TF 

91,60,10 
X2f97 

DO 

A 
A 

X2-7, IN7-2 
Ml, IN5 

S 

CM 

M1,0UT4 
Ml, 0,9 

BE 
TF 

NEXT7+36 
0UT11.ZER0-4 

MEXT7 

MM 
SF 

Ml, 24, 10 
96 

TF 
A 

TIFH-7,99 
TIFH.X2 

A 
S 

TIFH,DE3 
TIFH^TFH 

TF 
TF 

Ml,ZER0-7 
DE3,ZER0-1 

TF 
S 

X6, IN8 
X6,X5 

TF 
LD 


0UT9-2,X6 
_89*J(6_ 


D 
AM 


TF 
C 


93,TIFH 
_88^_ 


OUT8,87 
X6,ZER0-6 


BNN 
CHECK  CM 


BNE  CHECK 

CM   TIFH-7,1,8 


TESTK 
IN9.7.10 


BE 
B 


BREAK 
CONT 


TESTK  CM 
B£_ 


BREAK 


BNC 
BC1 


IN9,9,10 

CONT 

2F0UT 
RETURN 


B 
RETURNTFM 


SEARCH 
BRCH+6. SETUP 


TF   0UT10-2,X6 
B AROUND+1 2 


CONT 


C 
BE 


IN4,0UT3 
AROUND 


BNC 
BNC 


2CUDBE 
1CUDBE 


B 
CUDBE  TF 


AROUND 
DELTA-13, 


PMTF+11 


AROUNDA 
A 


BD   MOUT, ITTFC 

TF   DELTA-13, ORG+11 


0UT10-2,X6 
0UT11-2.X6 


LOAD 


A  0UT12-3,X6 

A  X5,X6 

BC1  SKIP 

TF  .I1.0UT8.2 


MM 
SF 


TIFH,6,10 
90 


"AM     94,5 

J^JOP   93^T1-10,7 


L0AD1    NOP   TI ,93, 2 

JLE NUMBER+11,1  nani+6,11 


AM      LOAD+6,10 

m i  pad  i  +^m 


AM      LOAD  1-1, 10 
TFM    LQADI+l^ft/in 


TFM   LOADI-11,21,10 
SKIP IE OUT!,  I  Nil 


TF 
J£_ 


TF 
JLF_ 


0UT2, IN2 
0UT3, IN4 


0UT4, IN5 
OUTS, IN6 
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AM      0UT17,1,8 

C         OUT9-2,ZERO-6 


BE 


*+36 
TIFH-2.5 


8 


A        TIME,TIFH-3 
JLE QUT7lz*,TIFH-7 


AM 


IN7,480,8 
HFRF 


CM     HR12+35,2,10 
BNE   *+36 


TF      TFH,X2 
JLE QUT6-4,X2-7 


SF      TIFH-6 

MM     TIFH-3,60,10 


AM     96,5,9 
J£ QUT7-2,95 


TF      TIFH,ZERO 

11 tifh-q,zero-8 


SF      X2-6 

MM      X2f6Qf1Q 


AM  93,5,9 
TF   QUT6-2fq2 


BRCH   BNC2SETUP,,2 
WRITE  WNCDOUT1-4 


TFM  ITTFC,9,10 
B   START 


END   TF  OUT19-2,ZERO-6 
, TFM  X5.Q, 8 


AM   I N7, 1200, 8 
3. HERE 


SM 
_B_ 


IN7,240,8 
HERE 


SEARCHCM 

m. 


OUTIO-2,25,8 
CLEAR 


CM   ITTFC,9,10 
BNE  CLEAR 


TF  0UT30,0UT1 
TF   OUT31,OUT2 


TF  0UT32,0UT3 
IE 0UT33.0UT4 


AM 
LQQP  TF 


TF  0UT34,0UT5 
ID 0UT40, CODES 


0UT42,1,10 
.RQA.ZERQ-3 


TF 
ACUM  A 


TF 
JIE- 


0UT21,YT 


ME 
YAMT 


0UT23-1, 
TD  0UT23,C0DES-2 
IE ni)T?5,QUT17 


AM  0UT25,1,8 
WNCD0UT19-13 


TF   TIME,ZER0-2 
IE YTIME.ZER0-2 


TF   YAMT,ZER0-4 

IE TFH.7FR0-4 

TFM  OUT42,0,10 
TFM  ITTFC.0.10 


B 

RFADERBD 


FIRST 
HR12, IN12-3 


BNF  *+36, IN7-3 
_££ LN7 


HERE 


AM      I N7, 1440,8 

SE IN7-3 


LD 
_DJ4_ 


92, IN7 

91.60,10 


TF 
_B_ 


X2,97 
JQD 


HR12      SF      IN7-3 

ID *+23,  INI  2-3 

CM     *+8, 01000, 710 
BNE  *+36 


£1 


AM 
_C_ 


TIA,ZER0-6 

RQA, I  2.7 

T I  A, 1,8 
TIA,T2,7 


BE 
_C_ 


_GIL 


BNE 
JLE_ 


NEW 
TIAr 


Nil 


ACUM 
RQAH.RQA 


TF   K2+11,ACUM+1 1 
IE C3+11,ACUM+11 


TF   SUB+11,F1+11 
TF   P3.RQA 


C3 


A 
A. 


REAL 


ROAH, 12,7 
B3,I2,7 


B3, I  2,7 

B3.R0A 


BNH 
J£_ 


SMALL  AM 
F2    C 


SMALL 
R0A,B3 


ANY, 1,8 
T2tANYt2 


BNE 
,AM_ 


BB 


AM 
AM 


BB 
REAL+11, 


10. 


F2+6,10 
TIA.1.8 


SUB 


C 
BNE 


TIA,T2,7 
CI .  ._2 


HI 


AM 
AM 


SUB+11,10 
K2+11.1Q 


AM  C3+11.10 
POT   C    SUB+11, LOAD  1+6 
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BE 
B 

XX 
C1..2 

C1 
C2 

C 
BE 

INT1+4,TIA,2 
RE 

NEW 

B 
AM 

C3 
ACUM+11,10 

AM 
C 

F1+1 1,10 
TIA. INT1 

BNE 
C 

ACUM 
F1+11.L0AD1+6 

BE 
B 

XX 

GO 

RE 

TFM 
TFM 

SUB+18,K2 
P0T+30.K2 

TF 
TF 

RESUB+11,SUB+1 1 
REK2+11,K2+11 

TF 

TF 

REC3+11,C3+11 
TIAH.TIA 

XX 

B 
LD 

K2 
98,R0A 

JJ 

D 
AM 

93,INT1,7 
95.5 

TF 
AM 

HOLD, 94 
CNT.24 

TF 
B 

*+1o,CNT 
0..2 

STORE 

NOP 
TF 

0,0 
OUT35.HOLD 

B 
TF 

BET 
OUT36.HOLD 

B 
TF 

BET 
OUT37.HOLD 

B 

TF 

BET 
OUT38.HOLD 

B 
TF 

BET 
OUT39.HOLD 

B 
TF 

BET 
0UT40-1.H0LD 

BEH 

B 

TFM 

BEH, ,2 
CNT.STORE-24 

AM 
TF 

0UT17, 1,8 
0UT41.0UT17 

WNCDOUT30-4 
TFM  BEH-6, FINAL 

BET 

TR 
AM 

OUT35-5,RECORD-5 
NUMBER+6,4 

NUMBERCM 
BH 

1 NT1 ,0, 2 

FINAL 

C    TIA,C1+6, 11 
BNE  RESET-12 


RESET 


TF  ROA,ROAH 

AM  JJ+11.4 

B  XX 

AM  C1+6.4 


TF 
TF 


ROA,ROAH 
TIA.TIAH 


REK2   TFM  K2+11,0 
RESUB  TFM  SUB+11.0 


REC3   TFM  C3+1 1.0,7 
TFM  SUB+18.C1 


IhM  P0T+30,C1 
TFM  REAL+11. 12 
TFM  F2+6,T2 
AM   JJ+11,4 


TFM  ANY, 0,8 

B    C3-12 

FINAL  TF 


TF 


R0AH,ZER0-4 
TIAH,ZER0-6 


TF   ANY,ZER0-6 
TFM  F1+11rT2 


TFM  ACUM+11, I  2 
TFM  REAL+11. I  2 


TFM 
TFM 


SUB+18,C1 
P0T+30,C1 


TFM  NUMBER+6, INT1 
TFM  JJ+1 1. INT1 


TFM  F2+6,T2 
TFM  C1+6, INT1+4 


TFM  BEH-6, BEH 
TFM  CNT.STORE-24 


AM  0UT17, 1 ,8 
TF   0UT41,0UT17 


TD   0UT40, CODES 
WNCD0UT30-4 


CLEAR  TR   OUT35-5,RECORD-5 
TF   I  1 PZEROr  2 


C   *-6,L0AD+6 
BE  *+36 


AM  *-30,10 
B   *-48 


TFM  *-54, I  1 
TFM 


L0AD1-1 , T1-10 


TFM  LOAD+6, I  1 
TFM  L0AD1+6,T1 
TFM  LOADI+1,41,10 
TFM  L0AD1-1 1 .41 . 10 


C    IN4,0UT3 
BE   RETURN 


CM   ITTFC,9,10 
BE   MOUT 
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SETUP 


B 

TFM 


RETURN 

99.60.10 

99,0UT6-2 


BH  HOUR 

TF   TT-2f0UT6-4 


TFM 
B 


TT,30,10 
*+36 


HOUR  TF  TT-2,0UT6-4 

TFM  TT,0,10 

TFM  BRCH+6,TRY 

TF  A-2.TT-2 


SETT1 


TF 
C 


BH 
BE 


A,TT 
Af0UT6-2 


FIXT1 
SETT2 


A 
B 


A, Til, 7 
SETT1 


FIXT1 


TF 
S 


99, A 
99*0UT6-2 


SETT2 


TF   TCK1.TI 1 ,27 

S    TCK1,99,2 


AM 
CM 


CNTB,1, 
CNTB^2A 


10 
10 


BE 
BH 


FIXT2 
Fl  XT4 


FIXT3  AM  FIXTI-13,2 

AM  Fl XT1+30, 2 

AM  FIXT1+35,2 

AM  FIXT1+42,2 


FIXT2 


B 
CM 


SETT1-24 
TI3. 60.10 


BNE  FIXT3 

TF   TCK3f0UT6-2 


FIXT4  TFM  CNTB,0,10  ■ 
TFM  FIXT1-13,TI1 

TFM  FIXT1+30,TCK1 
TFM  FIXT1+35.TI1 

TFM  FIXT1+42,TCK1 
TF  B.0UT6-2 

TF  B-2,0UT6-4 

CM  B-2f24f 10 


CK1 


CK2 


TF2 


CK3 


TF3 
"CRZT 


BRDG 
GOON 


BNE 
TFM 


WRITE 
B-2.0.10 


TRY 


B 
LD 


DM 
AM 


WRITE 
98r0UT8 


TRY1 


TF 
A 


93,60,10 
92+3- 


PM,96 
AMT1f IPM 


A 
A 


AM 
AM 


AMT2, IPM 
AMT3, IPM 


TCK1,1,10 
TCK2.1.10 


AM 
AM 


CM 
BNE 


TCK3,1,10 
Bf  1f  10 


B,60,10 
CK1 


AM 
TFM 


C 
BNE 


B-2,1,10 
B,0,10 


TCK1,TI1 
CK2 


TD   DIGIT, CODES 
TFM  TCK1,0,10 


SF 
AM 


AMT1-4 
AMT1.S 


TF1    TF   0UT13,AMT1-1 
TF   AMT1,ZER0-4 


C 
BNE 


TCK2,TI2 
CK3 


TD   DIGIT, CODES 
TFM  TCK2,0,10 
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APPENDIX    1    CONTINUED 


TD      D1GIT-1 ,ZERO 
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Read  header 
card 


Read  data 


Flag  the 
input  data 


Has  the  last 
card  been 
processed? 


G> 


PUNCH 


£} 


Is  the  Input 
precipitation 
depth  less  than 
the  previous 
depth? 


Set  previous 
depth,  X5=input 
depth 


11s  input  in  same 
month  as  previous 
input? 


y 


0    0-^0 


Is  this  the  first 
input  card? 


<D 


/Are  amounts  for 
regular  time  periods  and/or 
^maximum  intensity  parts 
being  used? 


<D 


Move  monthly 
summary  data 
to  the  output 
area 


Step  card 
number 


PUNCH 


.,_i 


Add  monthly 
totals  to 
yearly  totals 


Zero  monrbJjf 
and  daily 
totals 


Is  the  last  card 
the  end  of  a  storm? 


Are  amounts  for 
regular  time 
periods  being  used] 


Yes 


Are  maximum 
intensities  beinj 
\  determined? 
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(Is  this  the  last 
monthly  total  in 
the  year? 


fAre  amounts  for 
regular  time  periods 
and/or  maximum  inten- 
sity parts  being  used?/ 


0* — -©  ©     © 


DE3  =  2400  -  the 
previous  time  of 
day 


Zero  previous  time 
of  day 


Is  input  during  a 
leap  year? 


Transfer  yearly 
totals  to  output 
area 


Step  the  card 
number 


PUNCH 


Initialize  the 
program  for 
another  data 
group 


-O 


<D 


Set  days  in 
February  • 
28 


Set  days  in 
February  » 

22- 


Set  Ml  =  number  of  days  from 
first  of  previous  input  month 
to  first  of  present  input  month 


Is  the  input  from  a  chart 
reader? 


K^D 


Get  time  of 
day  in  hours 
and  fractions 
of  hours 


Get  input  time  as  number  of  minutes  from 
midnight  until  the  time  of  the  reading  by 
adding  to  the  chart  reader  time  the  correct 
number  according  to  the  digit  in  column  31 
of  the  input. 


Add  to  Ml  the  day 
of  input  and  sub- 
tract the  prior  day 
of  input 


Get  the  time  of  day  in  hours  and 
fractions  of  hours  by  dividing  by  60 
the  number  of  minutes  from  midnight 
until  the  reading  was  taken. 


Zero  daily 

precipitation 

total 

■ 

' 

Set  TIFH  ■= 
ML  multiplied 
by  24 

Time  interval,  TIFH  - 
TIFH  +  DE3  +  time  of  day 
-  previous  time  of  day 


*Zero  Ml  and  DE3 


Amount  of  precipitation  for 
time  interval  -  X6  -  input 
depth  -  X5 


Move  X6  value  to  output 


Compute  the  precipitation 
Intensity 
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(Does  X6  =  0? 


fboea   the  code 
Indicate  this  aB 
the  first  of 

^another  storm? 


(Is  the  input  in  \ 
the  same  month  as  1 
previous  input?    J 


Are  amounts  for 
regular  time 
periods  and/or  max- 
imum intensity  part^ 

being  used? 


(  V  /  Does  the  storm \ 

V     J\   cross  the  month?  ) 


Are  amounts  for  regular 
time  periods  being 
determined? 


Yes 


Are  maximum  intensities 
being  determined 


Add  X6  to  storm 
total 


X 


Add  X6  to  daily  and 
monthly  precipitation 
totals 


|  X5  -  X5  +  X6  | 


Set  storm  total  =  X6 


Are  maximum  inte 
^belng  determined 


gies    y_ ^Q 


Store  Intensity  for 
maximum  intensity 
search 


Compute  number  of 

minutes    in   time 
interval 


Sum  total  minutes  in 
storm  and  store  for 
maximum  intensity  search 


Set  maximum  intensity 
ending  for  test  in  search 
part  of  program 


X 


Step  these  storage 
locations 


Transfer  the  identification  and 
date  to  output 


|Step  the  card  number     | 


Bet  time  of  day  and 
time  Interval  in  hours 
and  minutes  and  transfer 
them  to  output  area 


X 


Are  amounts  for  regular 
time  periods  being 
determined? 


29 


G> 


Is  the  storm  accumulation 
less  than  0.25  inches? 


Is  this  the  first 
input  card? 


Transfer  the  storm 
ending  date  and 
identification  to 
ou t put  area 


Zero  intensity  and 
time  accumulations 


Add  first  storm 
intensity  to 
intensity  accumula- 
tion 


Add  1  minute  to  time 
accumulation 


Does  the  time 
accumulation  =  the 
accumulated  time  at  the 
end  of  the  storm  time 
interval? 


Does  the  time  accumulation  =  the 
first  maximum  intensity  time 
interval? 


Clear  the  maximum 
intensity  output 
area  


<D 


Set  to  store  values 
from  the  next  storm 


Is  the  input  in  the 
same  month  as  the 
previous  output? 


Is  this  the  first 
input  reading? 


o- 


Step  to  the  next  time 
interval  in  the  storm 


Does  the  time  accumulation  =  the 
first  maximum  intensity  time 
interval?  


L-0 


Does  the  first  maximum  intensity 
time  interval  ■  the  duration  of  the 
.storm? 


o 
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OH 


Set  B3  and  ROAH  =  intensity 
accumulation 


Add  proper  storm  intensity  to 
ROAH 


Add  proper  storm 
intensity  to  B3 


Subtract  from  B3 
the  first  added 
intensity 


'Does  T-i  =  the  time 
accumulated  at  the  end 
of  the  time  interval 
having  the  intensity 

.being  subtracted? 


Step  the  subtract 
instructions  and 
test  to  the  next 
time  interval  in 
the  storm 


Add  1  minute  to  the 
time  accumulation 


Divide  the  intensity 
accumulation  by  the 
number  of  minutes  in  the 
maximum  intensity 
interval 


<D 


Store  the  maximum  intensity  in 
the  proper  output  location 


T 


Is  this 
stored? 


the  sixth  intensity 


Is  it  the  twelfth 
intensity  stored? 


D 


Step  the  card 
number 


PUNCH 


Clear  the 
output  area 


Is  it  the  last  intensity  \ 
that  will  be  determined?  ) 


Make  setup  for 
next  maximum 
intensity 


Reset  all  values 
and  instructions  for 
the  next  storm 


Step  the  card 
number 


— * t_jT\ 

PUNCH    f— M      ) 


Does  the  time  accumulation  =  the 
time  accumulated  at  the  end  of  the 
storm  time  interval  from  which 
intensities  are  being  added? 


■*-?- 


Has  the  next  maximum 
intensity  time  interval 
been  previously  completed? 


Does  the  time  accumula- 
tion =  the  time  interval 
for  which  the  next  maxi- 
mum intensity  is  to  be 
determined? 


Step  the  add  instructions  and 
test  time  accumulation 
instructions  to  the  next  time 
interval  in  the  storm 


Was  the  last  time  interval  used  the 
last  one  in  the  storm? 


Store  values  for  starting  of 
next  maximum  intensity 
determination 
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Does  the  storm  start 
before  the  half  hour? 


Set  test  time 
on  hour  before 

storm  started 


-Hset 


Set  test  time  on 
half  hour  before 
storm  started 


A  =  test  time 


Compare  A  with  the  time  of 
day 


Subtract  from  A  the  time 
of  day  plus  the  regular 
time  period.   This  quantit  i 
is  the  period  time  check. 


Step  to  next 
regular  time 
period 


Does  the  third 
regular  time 
period  =  60? 


Set  time  check  for 
this  period  equal  to 
time  of  day 


I  Zero  CNTB,  reset  time  check 
to  first  regular  time  period 
and  transfer  time  of  day 
to  B. 
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0 


Get  precipitation  intensity  in 
inches  per  minute  -  IFM 


Add  IPM  to  each  of  the  three 
regular  period  amounts 


<D 


Add  1  minute  to  each  regular  time 
period  check  and  to  B 


Is  the  first  regular  time 
period  completed? 


1  Set  DIGIT  -9 


Zero  first  time 
period  check 


Transfer  amount  in 
period  to  output 
area 


Zero  amount  for  next 
accumulation 


Is  the  second  regular  time 
period  completed? 


I  Set  DIGIT  -  9   I 


Zero  second  time 
period  check 


Transfer  amount  in 
period  to  output 
area 


Zero  amount  for  next 

accumulation 


Is  the  third  regular  time 
period  completed? 


Set  DIGIT  -  9 

' 

Zero  third  time 
period  check 

f 

Transfer  amount  in 
period  to  output  area 

OH 


Zero  amount  for 
next  accumulation 
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»{  Does  DIGIT  =  9?   ) 


Put  time  of 
period(s) 
ending  in 
output  area 

I 

MHM 

V 

Set  the  code 
and  step  the 
card  number 

' 

f 

PUNCH 

' 

f 

Zero  output 
area 

■ 

t 

Set  for  next  output 
area  if  another 
period  can  be 
completed 

Does  B  =  time  of  day  of 
input  reading? 


Does  B  =  2400? 


Set  3  -  0000 


k. 


<D 


Have  any  time  periods  been 
completed  since  the  previous 
input  reading? 


Set  DIGIT 


Set  for  original 
output  area 


G 


^-O 
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O-H 


Do  all  regular  time  period 
accumulations  =  0? 


Yes 


Transfer  amounts  and 
B  to  output  area 


Zero  accumulations 
and  time  period 
checks 


PUNCH 


Initialize  the  amounts  for  regular 
time  periods  part  of  the  program  for 
the  next  event 


Is  the  input  in  the  same 
month  as  previous  input? 


Are  maximum  intens 
being  determined? 
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A 

V 
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/ 
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K: 
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